Anomalous Pinning Fields in Helical Magnets: Screening of the Quasiparticle 
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The spin-orbit interaction strength Qbo in helical magnets determines both the pitch wave number 
q and the critical field Hci where the helix aligns with an external magnetic field. Within a standard 
Landau-Ginzburg- Wilson (LGW) theory, a determination of Qao in MnSi and FeGe from these two 
observables yields values that differ by a factor of 20. This discrepancy is remedied by considering 
the fermionic theory underlying the LGW theory, and in particular the effects of screening on the 
effective electron-electron interaction that results from an exchange of helical fiuctuations. 
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Chiral itinerant ferromagnets such as MnSi [3, Q and 
FeGe Q have recently attracted considerable attention. 
Both of these systems crystallize in the cubic B20 struc- 
ture, which lacks inversion symmetry, and as a result 
spin-orbit coupling effects are important for the mag- 
netic properties. They both exhibit spiral or helimag- 
netic spin order at low temperatures (below about 28.5 
K at ambient pressure in MnSi, and below about 279 
K in FeGe, respectively), which is believed to be gener- 
ated by a Dzyaloshinski-Moriya term 0, 0] in the free 
energy. The pitch wavelength of the helix is large com- 
pared to a microscopic length scale; this reflects the weak- 
ness of the spin-orbit interaction. The pitch wave vector 
is q = 0.035 in MnSi and g w 0.009 l^^ in 
FeGe 3]. In other parts of the phase diagram in MnSi, 
striking non-Fermi-liquid behavior has been observed in 
low-temperature transport measurements . 



In an isotropic electron system there would be no pre- 
ferred direction for the pitch vector of the helix. In real 
materials, the underlying crystal lattice pins the helix. 
The terms in the free energy that cause this pinning are 
of higher order in the spin-orbit interaction, and hence 
represent an energy scale that is even weaker than those 
that lead to the formation of the helix. In MnSi, the pin- 
ning is in the (1, 1, l)-direction (or equivalent); in FeGe, 
it is in the (1, 0, 0)-direction (or equivalent) close to the 
transition, and in the (1, 1, l)-direction at lower temper- 
atures. An external magnetic field makes it energetically 
favorable for the helix to align with the field, and this 
competes with the crystal-field effects. As a result, upon 
applying a magnetic field in, say, the (0, 0, l)-direction 
to MnSi in the helical phase, the pitch vector rotates 
away from (1,1,1) until it aligns with the field direction 
at a critical field strength H = Hd- This field strength 
marks the boundary of the so-called conical phase, which 
is characterized by a homogeneous magnetization super- 
imposed on the helix that is aligned with the field. Upon 
further increasing the field, the homogenous component 



of the magnetization increases, and the amplitude of the 
helix continuously decreases until it disappears at a sec- 
ond critical field, i?c2, where the system enters a field- 
polarized phase with a homogeneous magnetization. 

The above considerations make it clear that one can 
obtain a measure of the spin-orbit interaction strength 
from measuring either the pitch wave number or the 
critical field Hd- It is a puzzling, but overlooked, fact 
that interpreting the results within the existing theoreti- 
cal framework yields values for the spin-orbit interaction 
strength that differ by a factor of about 25. It is the 
purpose of the present Letter to resolve this discrepancy. 

In order to frame our discussion of these various effects, 
let us consider the Landau-Ginzburg- Wilson (LGW) the- 
ory that has been commonly used to describe helical mag- 
nets Q. If the phase transition is either continuous or 
weakly first order, then the classical behavior of the sys- 
tem close to the transition can be described by an action 
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Here Sq is the usual action for a classical Heisenberg fer- 
romagnet 0], 



Sq 



dx 



(lb) 

Here M{x) is the three-component order parameter 
whose expectation value is proportional to the magneti- 
zation, ro is the bare distance from the critical point, and 
a and u are parameters that depend on the microscopic 
details of the system. (VTVf)^ is a shorthand notation for 
J2i j diMj d^M^ . So is invariant under separate rotations 
in order-parameter space and real space. 

Sc is the leading chiral term induced by the spin-orbit 
interaction 0,0], 
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J dx M{x) ■ (V X M{x)). 



(Ic) 
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The coupling constant c is proportional to the dimen- 
sionless spin-orbit interaction strength ^so, and on di- 
mensional grounds we have c = a/cpffso- Note that 
is still invariant under joint rotations in order-parameter 
space and real space, but not under spatial inversions. 
This terms can therefore be present only in systems that 
are not inversion invariant. The chiral nature of the curl 
produces the helical ground state, and the handedness of 
the helix depends on the sign of c. We assume c > 
without loss of generality. 

5'cf is the largest term that describes the crystal-field 
effects that couple the magnetization to the underlying 
lattice. For a cubic lattice, a representative contribution 
to Scf reads [sl 
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The coupling constant h is quadratic in g^o and given by 
h = a' gl^ with \a'\ w a. (Here and it what follows we ig- 
nore factors of 0(1).) 5cf breaks the rotational invariance 
and is responsible for pinning the helix. The direction of 
the pinning depends on the sign of b. 

Finally, we have a term that couples an external mag- 
netic field H to the magnetization: 



Sh^ j dx H{x) ■ M{x). 



(le) 



In the absence of both an external field and any cou- 
pling to the underlying lattice, it is easy to see that the 
Eqs. (jlbllcp lead to a helical ground state: 

M{x) = mi [ei cos(<7 ■ x) + €.2 sm{q ■ x)] . (2) 

Here the unit vectors ei, 62, and q — q/\q\ form a right- 
handed dreibein. The amplitude of the helix is given by 
TOi — y^—{t — aq'^)/u. The pitch vector q points in an 
arbitrary but fixed direction, and in a mean-field approx- 
imation its modulus is given by q — c/2a + 0{g^^). q is 
small compared to the Fermi wave number by virtue 
of the smallness of gso- In MnSi, kp « 3.6 |T3|, so 
g/fcp « 0.01. Assuming the same value for kp in FeGe, 
we have q/kp « 0.0025. The value of gso, which is equal 
to q/kp within a factor of 2 11], is thus 
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0.01 



(MnSi) 



0.0025 (FeGe). 



(3) 



As discussed above, a magnetic field tends to align the 
helix away from the pinning direction that is ultimately 
determined by the spin-orbit interaction, and hence the 
magnitude of the field necessary to depin the helix pro- 
vides another estimate for gso- In MnSi, the best stud- 
ied helimagnet, the helix in zero field is pinned in the 
(1,1, l)-direction, which implies that the coefficient b in 
Eq. (fid)) is negative (6 > leads to pinning in the (1,0, 0- 
direction) [8]. At ambient pressure, and not too close to 



the transition temperature, the experimental value for 
the field Hd defined above, where the pitch vector q 
aligns with the field direction, is Hd ~ O.IT [13]. In 
the same region, the experimental value for the field Hc2, 
where the helix vanishes, varies between 0.4 T and 0.55 T. 
Together with the corresponding experimental results for 
FeGe Q, we thus have a ratio 



Acxp — Hcl/Hc2 




(MnSi) 
(FeGe). 



(4) 



Comparing these experimental results with theoretical 
estimates leads to a puzzle. The Eqs. ([T|) can be ana- 
lyzed in detail to yield the critical fields Hd and Hc2 
HQ, but for our present purposes the following sim- 
ple considerations suffice. Hd is roughly determined by 
the magnetic energy, given by Sh, being equal to the 
pinning energy, which is given by Sd- Hc2 is roughly 
determined by the magnetic energy being equal to the 
chiral energy, which is given by Sc- For a homogeneous 
magnetic field, the coupling in Sh is to the homogeneous 
magnetization, uiq = J dx M{x) = xH , with x the ho- 
mogeneous magnetic susceptibility. The magnetic energy 
is thus of 0{H^), Sh = xH^- In Eq. ([ld|, the gradient 
squared is on the order of g^, and the magnetization is 
on the order of the amplitude of the helix, mi. For Hd 
we thus obtain the estimate 

Hd ~ 9so mi q^/a/x- (5a) 
Applying an analogous estimate to Eq. (fTc]). we obtain 

Hc2 ~ mi^ afcpgso q/x- (5b) 

All quantities whose estimates might be questionable 
thus drop out of the ratio A, and we have the theoretical 
result from the bare LGW theory 
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Comparing Eqs. dl]), (HI), and ©, we see that the experi- 
mental values for A are larger than the theoretical expec- 
tation by about a factor of 20 in MnSi and 40 in FeGe. We 
will now show how this discrepancy can be resolved by 
considering the screening of the effective electron-electron 
interaction that results from the exchange of helical fiuc- 
tuations. 

We first need to discuss the nature of the dominant 
fiuctuations in a helical magnet. The helical ground state 
represents a spontaneous breaking of translational invari- 
ance, and therefore leads to a Goldstone mode or heli- 
magnon. For the LGW action written above, and as a 
function of the wave vector A; for )fc) <C g, the frequency 
of the helimagnon reads [li[3 



wo(fc) 



akf.+2\b\kl/3 + akl/2q^. 



(7) 



Here k = (A;||,fe) has been decomposed into components 
parallel and perpendicular to the pitch wave vector q. 
We have assumed 6 < 0, as appropriate for MnSi [1], and 
we have neglected corrections of 0{b) = 0{gl^) to the 
coefficients of fc^ and k\. For & = 0, that is, if we neglect 
all effects of 0{g^^), the helimagnon frequency squared 
lacks a contribution proportional to A;^ . This is a result 
of rotational invariance. The crystal-field term 5'cf in 
the action breaks this invariance, which leads to a mode 
that is still soft (because the translational invariance is 
still broken), but has a fc^-term with a small prefactor. 
It is the generalization of the well-known magnons in 
ferromagnets and antiferromagnets that have a quadratic 
and linear dispersion relation, respectively. 

The spin model described by the Eqs. H]) can be under- 
stood and derived as an effective theory that results from 
an underlying fermionic action. The technical procedure 
is to single out the magnetization by either performing 
a Hubbard- Stratonovich decoupling of the spin-triplet 
interaction, or by constraining the appropriate combi- 
nation of fermion fields to an auxiliary composite field 
whose expectation value is the magnetization 
Integrating out the fermions then yields the spin model, 
with the coefficients of the LGW theory given in terms of 
localized fermionic correlation functions. Conversely, in- 
tegrating out the magnetization yields an effective theory 
of electronic quasi-particles that interact via an exchange 
of helimagnons. This effective interaction was derived 
and discussed in Ref. [3. For small wave numbers, the 
leading contribution to the effective potential is 

V{k,in;pi,p2) = Vaxik,ifl)-/ik,pi)-/{-k,p2). (8a) 

Here Vq — A^q^/8(m*)^ with A the Stoner gap (i.e., the 
splitting of the two electron bands that results from the 
magnetization in a mean-field approximation) and m* 
the electron effective mass, denotes a bosonic Matsub- 
ara frequency, and 



xik,in) = 



1 



1 



2A^F 3fc| uj^ik) - {iiiy^ 



(8b) 



with Np the density of states at the Fermi surface is the 
helimagnon susceptibility. The leading contribution to 
the vertex function 7 is given by 



7(fc,p) = iy{ki_ ■pi_)p\\/X. 



(8c) 



Here v is a dimensionless parameter that describes the 
coupling of the electrons to the lattice; generically, one 
expects V = 0(1) [3]. The effective potential is graph- 
ically depicted in Fig. [TJ notice that it depends on the 
momenta of the quasiparticles involved in addition to the 
transferred momentum. 

Due to the singular nature of the helimagnon suscepti- 
bility the effective interaction is long-ranged, and screen- 
ing has a qualitative effect. The leading effect of screen- 
ing is captured by the usual random-phase approximation 



k,iQ 



< 



J 



P2 -J:, ico^-in 



FIG. 1; The effective quasiparticle interaction due to heli- 
magnons. Note that the vertices depend on the quasiparticle 
momenta in addition to the helimagnon momentum. 



(RPA) 20]. The screened interaction is shown in Fig. 
[21 If one takes the resulting fermionic theory of quasi- 
particles interacting via the screened effective interaction, 
re-introduces the magnetization, and integrates out the 
fermions, one can study the effect of the screening on the 
LGW spin model. For our present purposes, the most 
important effect is a renormalization of the coupling con- 
stant in the crystal-field term Sd, Eq. (|ld[) . We find 



(9) 



While the value of the Stoner gap in MnSi is not well 
known, it is clear from the small value of the critical tem- 
perature that ep/A is large compared to unity. An band 
structure calculation in Ref. I21J yielded A w 3,000K, 
while the Fermi temperature is Tp « 150, 000 K ^10., ,22|. 
for a ratio ep/A w 50. The effect of the screening is thus 
large, and the sign of the effect is of fundamental im- 
portance. If 6 < 0, then the bare theory yields a helix 
pinned in the (f , f , f )-direction and the renormalization 
greatly enhances the coefficient and hence the pinning 
strength. If 6 > 0, then the bare theory predicts pinning 
in the (1,0, 0)-direction and the renormalization changes 
this to produce pinning in the (I, I, f )-direction. This 
is analogous to the effect of a strong electron-phonon 
coupling that can trigger a structural phase transition. 
The conclusion that the pinning will always be in the 
(1,1, f )-direction holds for systems where the transition 
is continuous or weakly first order. For strongly first or- 
der transitions the LGW theory is no longer controlled, 
and a gradient-free cubic anisotropy in the action (i.e. a 
term proportional to M^, which we have neglected in 
Eqs. IJ)) can invalidate this conclusion. These observa- 
tions are consistent with experimental results. In MnSi, 
the transition is continuous or very weakly first order, 
and the pinning is in (1, 1, f )-direction everywhere in the 
ordered phase [4]. In FeGe, the transition is strongly first 
order, and the pinning is in (1, 0, 0)-direction close to the 
transition, but switches to the (1, 1, f ) direction at lower 
temperatures 

From the considerations leading to Eqs. ^ we see that 
Hci is proportional to Vb while Hc2 is independent of b. 
The renormalized theory thus yields the following result 
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Pi+k P2-k p, + k p^-k p^ + k p^-ic p^-k 

FIG. 2: Screening of the effective quasiparticle interaction. 



for the ratio A = Hci/Hc2, 

Athco ~ gsoi^ep/A, (10) 

which replaces Eq. ([6]). With the numbers quoted above, 
and assuming « 1, we obtain values for MnSi and FeGe 
that are in agreement with the experimental ones given 
in Eq. (j4]) within a factor of 2. 

As a check, we finally discuss the absolute values of 
Hci and Hc2- Taking into account the renormalization 
of 6, Eq. ([9|), Hci is given by 

Hci ~ 9so miqv (ep/A) a/ x- (H) 

In the ordered phase of MnSi at ambient pressure, the 
susceptibility is observed to be roughly x ~ S^I/fceTc 
"M.- This is consistent with theoretical considerations 
231 ■ In a fully renormalized spin model, the gradient 
squared term in the action taken at the Fermi length 
scale must be roughly equal to the critical temperature, 
ak^m\/2 « T^. With these estimates, we obtain from 
Eq. (fTTj) Hci ~ 0.25 T, which is the correct order of mag- 
nitude. Similarly, from Eq. ()5bp we obtain Hc2 ~ 0.5 T, 
in agreement with the experimental value. 

In summary, we have pointed out that the standard 
LGW theory for helical magnets leads to a large discrep- 
ancy between the strength of the spin-orbit interaction in 
helical magnets as determined from the pitch wave num- 
ber versus the ratio Hci/Hc2 of the two critical fields. 
We have shown that a renormalization of the theory that 
results from the screening of the effective quasiparticle 
interaction resolves this puzzle and also leads to absolute 
values of the critical fields that are in good agreement 
with the experimentally observed values. 
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